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Abstract: Measurements of heavy flavor quark (charm and bottom) correlations in heavy ion 
collisions are instrumental to understand the flavor dependence of energy loss mechanisms in hot 
and dense QCD media. Experimental measurements of these correlations in baseline p+p collisions 
are crucial to understand the contributions of perturbative and non-perturbative QCD processes 
to the correlation functions and further help in interpreting correlation measurements in heavy ion 
collisions. In this paper, we investigate D-D meson correlations and D with one particle from D 
meson decay daughter correlations using PYTHIA Event Generator in p + p collisions at >/s = 200, 

500 and 5500 GeV. Charm/bottom events are found to contribute mainly to the away side/near 
side pattern of D-electron correlations, respectively. In the energy region of RHIC, D-D correlations 
inherit initial c-c correlations and B -* DX decay contribution is insignificant. Furthermore, Bottom 
quark correlations are suggested to be applicable at LHC energy, as the bottom contributions on D 
related correlations are relatively large. 

PACS numbers: 25.75.-q, 25.75.Cj 


I. INTRODUCTION 

Experimental observables indicate that strongly cou¬ 
pled Quark-Gluon Plasma (sQGP) is formed in the high 
energy heavy ion collisions W- Heavy quarks (c and b) 
serve as good probes to study the properties of sQGP, 
as charm and bottom can only be pair-produced in 
early stage hard scatterings due to the large mass (> 1 
GeV/c^) 0. Experimental measurements have shown 
that heavy quarks suppression at relatively high-pp is 
as large as that of light quarks 013 , which was incon¬ 
sistent with earlier theoretical expectation of the flavor 
dependence of parton energy loss. Thus it is crucial to 
understand the charm-medium interaction mechanism. 
Recently, the azimuthal correlations of heavy quarks are 
found to have the potential for distinguishing different 
energy loss mechanisms inside the hot medium ISE]. 
The theoretical prediction indicates that pure radiative 
energy loss does not change the initial angular correlation 
function in a significant way, whereas pure collisional en¬ 
ergy loss is more efficient at diluting initial back-to-back 
charm pair correlation, this could even lead to a peak in 
the near side at low-pp m- To approach the azimuthal 
correlations of heavy quarks in experiment, measurement 
of D-D correlations is the ideal candidate uni HI], as D- 
D correlations inherit most of charm pair correlations. 
However, due to the relatively small charm cross section, 
small hadronic decay branching ratio which is used to 
reconstruct D meson and limited signal to background 
ratio, the direct measurement of D-D correlations would 
be very challenging. One has to consider use D-X or X- 
D correlation measurement instead of direct D-D, where 
X is a decay daughter of D meson. 

Gonsidering only leading order perturbative Quantum 


Ghromodynamics (pQGD) processes, the initial charm 
pairs exhibit exactly back-to-back correlation. High or¬ 
der pQCD processes broaden the initial charm corre¬ 
lations. Gharm quarks finally fragment into charmed 
hadrons, the non-perturbative fragmentation processes 
further broaden the correlations. In the case of D-X or 
X-D correlation, decay kinematics and random combina¬ 
torial pairs would dilute or even destroy the correlation 
signal originating from initial charm pair correlations. It 
is crucial to understand these effects on the correlation 
function firstly in p -|- p collisions which serve as the 
baseline for heavy ion collisions. Furthermore, an exper¬ 
imental measurement of D-D and D-X {X-D) correla¬ 
tions in p -k p collisions would constrain the contributions 
from high order processes of pQGD calculations and non- 
perturbative fragmentation processes. In this paper, we 
utilize PYTHIA Event Generator (EG) to study D-D 
and D-X {X-D) correlations in p -k p collisions at Rel¬ 
ativistic Heavy Ion Collider (RHIC) and Large Hadron 
Collider (LHC) energies, ^/s = 200, 500 and 5500 GeV. 

The rest of the article is organized in the following way: 
Section H summarizes the details of tuning PYTHIA pa¬ 
rameters. In Sec. HI, we present the correlations results 
of D-D and D-X {X-D). The comparsion is made be- 
ween RHIC and LHC energies. Last a summary is pre¬ 
sented in Sec. IV. 


II. PYTHIA TUNE 

PYTHIAS (versions. 168) has been used in this study. 
In the following, we will use PYTHIA representing 
PYTHIAS. We tuned PYTHIA parameters to match 
the cc production cross-section as inferred from measure- 
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FIG. 1: (Color online) cc production cross section as a function of transverse momentum in p + p collisions at = 200, 500 
and 7000 GeV from STAR and ALICE experimental data (solid squares) and tuned PYTHIA (red lines). 


merit of and D* production m p + p collisions at 
y/s = 200 and 500 GeV measured in the STAR exper¬ 
iment at RHIC [nng and 7000 GeV measured in the 
ALIGE experiment at LHC El- The choice of modify¬ 
ing the strong interaction coupling constant (a^) value 
of final parton shower (TimeShower:alphaSvalue) and 
minimum invariant transverse momentum {pr) thresh¬ 
old for hard QCD process (PhaseSpace:pTHatMin) gives 
the best x^/ndf to the data. Figure 1 shows the cc pro¬ 
duction cross section as a function of px from STAR mea¬ 
surement and PYTHIA. In the PYTHIA calculation, all 
ground state charm hadrons (ZI°, D^, , D~, 

A); and A“) were added together in the rapidity window 
\y\ < 1 to obtain charm cross-section. The parameters of 
‘TimeShower: alphas value’ and ‘PhaseSpace:pTHatMin’ 
were set to 0.18 and 1.3 GeV/c for = 200 GeV col¬ 
lisions. These two parameters were set to 0.15 and 1.5 
GeV/c for ^/s = 500 GeV collisions; 0.15 and 2.8 GeV/c 
for 7000 GeV. The y^/ndf values are 1.31 for 200 GeV 
data, 2.13 for 500 GeV data and 1.95 for 7000 GeV data, 
respectively. As lack oi p p experimental data in top 
energy of heavy ion collisions at LHC (5500 GeV), the 
same setup of parameters as 7000 GeV was applied in 
PYTHIA simulation. 


III. CORRELATIONS 

D-D correlations are the ideal experimental way to 
study the azimuthal correlations of initial c-c pairs, as 
the charm quark mostly fragment to D mesons. In this 
paper, we use D to represent the combination of 11°, D^ ^ 
D^ and D to represent D~, D^. In experiment, D 
mesons can be reconstructed by the hadronic decay chan¬ 
nels. Figure 2 shows the D-D correlations vap p colli¬ 
sions at \/s = 200 GeV for (a) real pairs and (b) all pairs. 
Where real pairs mean D-D from c-c pairs, all pairs mean 


all possible D-D pairs in an event. A phase space cut of 
pseudorapidity {\ri\ < 1) was applied to match the typi¬ 
cal STAR experimental acceptance. The integral of all 
the correlation functions were normalized to one, as we 
are interested in the correlation patterns here. The D-D 
azimuthal correlation function exhibits a clear away side 
correlation inherited from the initial c-c pair production. 
In case we only consider the leading order (LO) pQCD 
contribution, the c-c would show a delta function like 
away side correlations. The next-to-leading (NLO) or¬ 
der pQCD contributions including radiative corrections, 
flavor excitation and gluon splitting broadens the corre¬ 
lation function. There are no exact NLO pQCD calcula¬ 
tions in PYTHIA, but the initial and final parton shower 
procedure emulates these effects. Thus some aspects of 
the multiple parton emission phenomenon could be well 
reproduced [m ng. The non-perturbative fragmenta¬ 
tion process further broadens the correlation function. 
In other words, the correlation function of D-D inher¬ 
its the initial c-c back-to-back correlation broadened by 
higher order pQCD and fragmentation processes. The 
results of all pairs in panel (b) are consistent with real 
pairs in panel (a), because most events only contain one 
D-D pair with phase space cuts ((|? 7 | < 1)) in p -|- p col¬ 
lision at RHIC energies. The experimental measurement 
of D-D correlations in p -|- p collisions is crucial to con¬ 
strain the higher order pQCD and fragmentation calcula¬ 
tions. Further in the heavy ion collision system, Au -|- Au 
for example, the initial away side correlation is expected 
to by the charm-medium interactions. Different energy 
loss mechanisms, collision energy loss or radiative energy 
loss, show dramatically different modification [8l|9]. The 
experimental measurement of D-D correlation in heavy 
ion collision system will help us to disentangle the ques¬ 
tion on charm quark energy loss mechanism. However, 
it is a quite challenging analysis, the statistics are lim¬ 
ited by the small c-c production cross section, the usable 
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FIG. 2: (Color online) The D-D, D-K, K-D, D-e and e-D correlations as a fnnction of relative azimuth angle A<() in p + p 
collisions at ^/i = 200 GeV in PYTHIA for (a) real pairs (pairs are from c-c pairs), b) all pairs (all possible pairs in an event). 
The phase space cnt is pseudorapidity (|p| < 1). An additional pr cut is applied to trigger particles. The integrals of correlation 
functions are normalized to one. 


hadronic decay branch ratio (e. g. K~ + tt^) and 

the limited signal over background ratio of reconstructed 
D mesons |12) . 


A. D-X (X-D) correlations 

Instead of direct D-D analysis, D-X or X-D correla¬ 
tions are other options to approach charm quarks cor¬ 
relations, where D-X means correlations of D{D) me¬ 
son and a decay daughter of D{D) meson (electrons, 
kaons, pions and so on). (First letter indicates the trig¬ 
ger particle.) The typical choices would be charged kaons 
and electrons. Charged pions are not good candidates, 
because of the copious production and resonance de¬ 
cay contributions to pions. Electrons could be divided 
into photonic and non-photonic electrons. Photonic elec¬ 
trons include those from 7 conversion and Dalitz decay, 
while non-photonic electrons include those from charmed 
and bottomed hadrons semileptonic decays [SI IZ]- Ex¬ 
perimentally, non-photonic electrons can be statistically 
subtracted from inclusive electrons li m. Thus in the 
PYTHIA simulation results in Fig. 2, we use non¬ 
photonic electrons for D-e and e-D correlations. Also 
shown in Fig. 2 are D-K and K-D correlations. Panel 
(a) shows the correlation function inheriting from initial 
c-c pairs. The away side correlations can be seen for 
all combinations. For D-D, D-K and D-e, the different 
widths of away side peak are due to the decay smearing 
oi D ^ K or D ^ e. The different widths of away side 
peak of D-D, K-D and e-D are because minimum trig¬ 
ger pt cut on D decayed particles (> 2 GeV/c) actually 
requires higher px D mother trigger particles. Panel (b) 


shows the correlation function for all possible pairs in an 
event. The large near side peak for D-K and K-D cor¬ 
relation are mainly due to the jet correlation, as most 
kaons are not from D decays. It suggests D - inclusive- 
hadron correlations are not good candidates to study the 
charm correlations. e-D correlation also shows a near 
side correlation pattern, which is due to the semileptonic 
decay of B contribution. The STAR measurement indi¬ 
cates that ratio of non-photonic electrons from B to those 
from D is 0.2 at pt = 2 GeV/c [T7] . The random pairs 
which do not originate from the same c-c pair cause addi¬ 
tional combinatorial background. It is partly the reason 
why the away side correlation width of D-e is wider than 
that of D-D. We also found tighter or looser cut on mini¬ 
mum transverse momentum of associate particle does not 
change the correlation patterns dramatically. 


B. Components of D-D correlations 

The D{D) mesons consist of directly produced D{D) 
and D{D) from B{B) meson decays. To study the B 
contribution to the D-D correlation function, we sepa¬ 
rate the D-D pair combination into four cases: 1) direct 
D and D, 2) D from B decays and D from B decays, 
3) D from B decays and direct D, 4) direct D, and D 
from B decays. The combination of case 2) -|- 3) -1-4) 
is generally called D-D pairs from B decays. Figure 3 
shows the decomposed D-D correlations in p -\- p colli¬ 
sions at \/s = 200 GeV ((a) and (b)) and 500 GeV ((c) 
and (d)) for two sets of trigger particle minimum pp cut. 
It can be observed that the B decay contribution to D-D 
depends on trigger particle pp cut and collision energy. 
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FIG. 3: (Color online) The D-D correlations from different sources as a function of relative azimuth angle Acp in p + p collisions 
at \/s = 200 GeV (panel (a) and (b)) and 500 GeV (panel (c) and (d)). The minimum trigger pr cuts are set to 2 GeV/c 
(panel (a) and (c)) and 3 GeV/c (panel (b) and (d)). The correlation functions are normalized by the number of D mesons. 


The higher trigger particle pt cut or collisions energy, 
the larger contribution from B decay. In comparison to 
direct D-D pairs, D-D pairs from B decays are insignif¬ 
icant. In p -I- p collisions at ^/s = 500 GeV, with p^*® > 
3 GeV/c cut, the D-D pairs from B decays are ~ 4% of 
the total D-D pairs. As shown in the panel (d) of fig. 3, 
the B contributions cause the small near side correlation 
pattern of the total D-D correlation. In p -h p collisions 
at s/s = 200 GeV, even with p^'® > 3 GeV/c cut, the 
contribution from B decays is still not sizable. It is also 
found that tighter or looser cut on minimum transverse 
momentum of associate particle does not change the con¬ 
clusion. It indicates the D-D correlations mostly inherit 
from c-c pair production at RHIG energies. The effect of 
B decays on D-D correlations is insignificant. 


C. D-e correlations 

H-electron {D-e) correlations contain contribution 
from different sources, thus it is more complicated than 
D-D correlation. The inclusive electrons consist of pho¬ 
tonic electrons and non-photonic electrons. The pho¬ 
tonic electrons are from 7 conversion (in experiment) 
and Dalitz decay (e. g. tt'^, p), where the non-photonic 
electrons are from semileptonic decay of charmed and 
bottomed hadrons. The non-photonic electrons can be 
statistically separated from the inclusive electron sam¬ 
ples in experiment laizi which makes the D and non¬ 
photonic electron correlation measurement possible. Fig¬ 
ure 4 shows the D^-e correlations from different electron 
sources in PYTHIA. Following the typical experimental 
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FIG. 4: (Color online) The D^-e~{D^-e^) and D^-e^{D^-e~) correlations as a function of relative azimuth angle At/i in p + 
p collisions at s/s = 200 GeV. Different lines represent differnt sources of electron sample. The minimum trigger pr cut is 2 
GeV/c and the minimum associate pt cut is 1 GeV/c for panel (a) and 1.4 GeV/c for panel (b), respectively. 


way, we chose from golden hadronic decay channel 
+ TT~. Based on the sources of electrons, we plot 
the the correlation of D°: 1 ) to electrons from charmed 
hadrons, 2) to electrons from bottomed hadrons, 3) to 
photonic electrons. Note that charmed hadrons contain 
the direct charm and decay contributions from bottomed 
hadrons, to be consistent with the experimental capa¬ 
bility. The inclusive electrons are simply the sum of all 
three cases mentioned above. As the low pp electrons 
are mainly from photonic sample, a minimum pp cut is 
usually applied on electrons. Panel (a) and (b) show the 
minimum pp cut of electrons = 1 and 1.4 GeV/c respec¬ 


tively. The initial c-c correlations contribute to D^-e~ 
and Z?b_g+ correlations which are shown in the left pan¬ 
els, where and D^-e~ can be interpreted as other 

contributions which are shown in the right panels. It is 
clear that the near side correlation pattern of D^-e~ is 
mainly from bottom events and the away side correla¬ 
tion pattern of D^-e~ is mainly from charm events. As 
discussed in reference [18], the azimuthal correlation of 
and non-photonic electrons allows the separation of 
charm and bottom production on a statistical basis. The 
- photonic e~ and - photonic correlation pattern 
are symmetric, because the Daliz decay or the 7 conver- 





















6 



A(t) (rad.) 


FIG. 5: (Color online) The D-D, D-e and e-D correlations as a function of relative azimuth angle A(j> in p p collisions 
calculated by PYTHIA at (a) >/s = 200 GeV and (b) 5500 GeV. The integral of correlation function are normalized to one. 
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FIG. 6: The ratio of direct D-D pairs to total D-D pairs as a function of trigger transverse momentum in p + p collisions at 
,/s = 200 and 5500 GeV. 


sion contributes equally to and e~. As the minimum 
Pt cut on electrons increases, the bottom contributions 
become larger. This can be explained by the increase of 
B ^ e to D ^ e ratio as a function of pp m The results 
of p + p at y/s = 500 GeV is similar to the 200 GeV re¬ 
sults. To investigate the c-c correlations by measuring D- 
e correlation, one has to stick to away side region which 
is dominant by charm contributions. The region of 2 


- TT would be fine inferring from PYTHIA simulation. 


D. RHIC versus LHC energy 

From RHIG to LHG, as the collision energy increases, 
higher order pQGD and bottom contributions to the cor¬ 
relation functions become more and more significant. In 
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Fig. 5, we compare the D-D, D-e (non-photonic elec¬ 
trons) and e-D correlations in p p collisions at s/s = 
200 and 5500 GeV. We chose these two collisions ener¬ 
gies, as they are the top heavy ion collision energies at 
RHIC and LHC, respectively. One can observe the near 
side correlations show almost the same magnitude as the 
away side correlations for D-D correlation function at the 
LHC energy. Furthermore, the higher order pQCD pro¬ 
cesses smear the away side correlations significantly: in 
the D-e case, the correlation function is almost flat in p -|- 
p collisions at ^/s = 5500 GeV. The large near side peak 
of e-D correlations are from B to e contributions. Simi¬ 
lar as Sec. HI B, we separate the D-D correlations into 
direct D and B to D correlations. It is found the near 
side correlations are from B to D contributions. Quanti¬ 
tatively, in Fig. 6, we show the ratio of direct D-D pairs 
to total D-D for 200 and 5500 GeV. It is observed that 
the D-D pairs from B contributions are relatively larger 
in LHC energy than RHIC energy. Up to trigger pT cut 
at 6 GeV/c, the D-D pairs from B contributions are less 
than 5% in 200 GeV collisions where the B contribution 
is ~I5% in 5500 GeV collisions. The PYTHIA results 
suggest that charm correlations might not be the best 
choice at the LHC energy, as the baseline measurements 
inp -\- p are significantly affected by B contributions and 
higher order pQCD. On the other hand, bottom corre¬ 
lations, especially B-B correlations are clean probes for 
both p -\- p and heavy ion collisions and possible to be 
done at LHC energies, due to the larger production cross 
section. 


IV. SUMMARY 

With PYTHIA (version 8.168) framework, we studied 
D-D and D-X (X-D) correlations in p -|- p collisions at 


^/s = 200 and 500 GeV. H-hadron correlation is found 
to be effected significantly by the jet correlations. D- 
D correlation is the ideal choice to approach initial c-c 
pair correlations. The bottom decay contributions to D- 
D pairs are less than 4% at p -I- p collisions in = 
500 GeV with trigger px cut at 3 GeV/c, thus it is not 
significant at RHIC energies. D - non-photonic electron 
correlation contains both contributions from charm and 
bottom events. The charm/bottom events dominate the 
near/away side correlation pattern respectively. D - non¬ 
photonic electron correlations can be used to determine 
the branching ratios for the charm and bottom decays to 
electrons. One can also study c-c correlations from the 
away side pattern of D-e correlation with a proper 
cut. The charm quark correlations are found to be less 
attractive at LHC energies, due to the relatively larger 
bottom and higher order pQCD contributions. Alterna¬ 
tively, the bottom quark correlations would be applicable 
at LHC. Experiment measurements on heavy quark (c 
and b) correlations are crucial to understand the charm 
correlations in p -|- p collision system, thus offering pre¬ 
cious data to constrain the pQCD calculations. Further¬ 
more, in heavy ion collisions, these measurements are 
important to investigate the charm-medium interactions 
and disentangle the contributions of different charm en¬ 
ergy loss mechanisms in the hot and dense medium. 
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